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We study the dynamical properties of the fully frustrated Ising model. Due to the absence 
of disorder the model, contrary to spin glass, does not exhibit any Griffiths phase, which has been 
associated to non-exponential relaxation dynamics. Nevertheless we find numerically that the model 
exhibits a stretched exponential behavior below a temperature T p corresponding to the percolation 
transition of the Kasteleyn-Fortuin clusters. We have also found that the critical behavior of this 
clusters for a fully frustrated q-state spin model at the percolation threshold is strongly affected 
by frustration. In fact while in absence of frustration the q — 1 limit gives random percolation, in 
presence of frustration the critical behavior is in the same universality class of the ferromagnetic 
q — I/2-state Potts model. 



o 

I. INTRODUCTION 

Glassy systems at low temperatures undergo a transition characterized by the freezing of structural relaxation, in 
which the system is trapped in a disordered metastable configuration. Already at temperatures higher than the ideal 
glass transition temperature To, a number of dynamical anomalies is observed. One of these anomalies concerns the 
relaxation functions of the system, which at high temperatures are characterized by a single exponential. Below some 
temperature T* higher than To, the long time regime of correlation functions, called a-relaxation, is well approximated 
by a Kohlrausch- Williams- Watts function, also known as "stretched exponential" , 

o 
o 



f(t) = foexp-(t/Tf (1) 



This behavior has been observed in many real glasses, such as ionic conductors, supercooled liquids and polymers 

A similar behavior has been observed in canonical metallic and insulating spin glasses, investigated by neutron 
and hyperfine techniques |8|-[T3|. These systems can be described by an Ising model, in which ferromagnetic and 
£SJ ', antiferromagnetic interactions are distributed in a disordered way on the edges of the lattice. The Ising spin glass 
■ undergoes a transition at some temperature T sg , called spin glass transition, analogous to the freezing transition of real 
glasses. Moreover, like in glass-forming systems, one observes a temperature value T* > T sg where non-exponential 
relaxation functions appear. This has been observed for the Ising spin glass in 2D by McMillan j^J, and in 3D by 
Ogielski 0. 

Several mechanisms have been proposed to explain the onset of non exponential relaxation functions like (|l]) in spin 
glasses, when the system approaches the glass transition from above. Randeria et al. [ [l6[ suggest that the temperature 
T* coincides with T c , the critical temperature of the ferromagnetic model. They base their conjecture on the presence, 
i-^j ■ in the spin glass, of non frustrated ferromagnetic-type clusters of interactions, the same that are responsible for the 
Griffiths singularity |l^] . The presence of non exponential relaxation in this approach is therefore a direct consequence 
O ■ of the quenched disorder. 

' On the other hand, it has been suggested [fj~8|— p0|| that in the spin glass the onset of stretched exponential relaxation 
^ , functions may coincide with the percolation temperature T p of the Kasteleyn-Fortuin and Coniglio-Klein clusters 
pl| , p2[ . These clusters can be obtained introducing bond with probability ps = 1 — e -2j/fc B T between nearest- 
neighbor pair of spins satisfying the interaction. In the spin glass the percolation threshold of these clusters is higher 
than T sg §§. Mor eover Stauffer j24j] observes non exponential relaxation in the d — 2 ferromagnetic Ising model, 
simulated by conventional spin flip, for temperatures lower than the critical temperature T c , which coincides with the 
percolation temperature T p . In d = 3 this behavior disappears, and the relaxation is purely exponential for all the 
temperatures. 

To gain more insight into the mechanisms which lead to the appearing of anomalies in the dynamical behavior, in 
this paper we study the fully frustrated Ising model p5|-p7| on a bidimensional square lattice, using the conventional 
spin flip techniques. In this model ferromagnetic and antiferromagnetic interactions are distributed in a regular way 
on the lattice (see Fig. [|), so that no unfrustrated cluster of interactions exists. Therefore in this model the appearing 
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of stretched exponential in the relaxation functions cannot be due to the mechanism of Randeria et al. We in fact 
find that the onset of stretched exponentials occurs at the percolation temperature T p , as we show in Sect. [n| 

The fully frustrated Ising model can be mapped onto a fully frustrated g-bond percolation by applying the Kasteleyn- 
Fortuin and Coniglio-Klein cluster formalism, where q = 2 is the multiplicity of the spins (see Sect. III). This model, 
which can be generalized to any value of q, is suitable to describe systems with geometrical frustration, and may give 
insight into the origin of the long relaxation decay in glasses, characterized by stretched exponentials |p8| . We study 
the dynamics of this geometrical model for q — 1 and q = 2, using the "bond flip" dynamics, that will be described 
in Sect. IV. We find again that the onset of stretched exponentials starts at the percolation transition T p (Sect. |v| ). 
Moreover we also find that the percolation transition is in the same universality class of the g/2-state ferromagnetic 
Potts model, as we show in Sect. [v|, in agreement with the expectation based on renormalization group |29|j . 

These results thus suggest that the presence of a percolation type transition may be responsible for the appearing 
of the "large scale" effects of frustration, among which there is the onset of various dynamical anomalies, such as 
stretched exponential relaxation functions |5(],^0| . 

In the frustrated g-bond percolation, the frustration is present at all length scale. To probe the effect of f rustr ation 
we have modified the model in such a way that only loops of length four are considered frustrated (Sect. VII). We 
find that the model for q = 1 has the same critical exponents of the ferromagnetic Potts model with spin multiplicity 
q = 1 (the random bond percolation model). Namely this local frustration does not change the critical behavior 
compared with the unfrustrated model. At the same time, although the dynamics of the model is influenced by the 
frustration constraint, in the long time regime the relaxation is purely exponential. 



II. THE RELAXATION FUNCTIONS OF THE FULLY FRUSTRATED ISING SPIN MODEL 

We simulate by conventional spin flip the fully frustrated Ising spin model, defined by the Hamiltonian 

n = -jJ2(^s i s j -i), (2) 

<«) 

where £y are quenched variables which assume the values ±1. The ferromagnetic and antiferromagnetic interactions 
are distributed in a regular way on the lattice (see Fig. |l|). 

We calculate the relaxation functions of the energy. Averages were made over 32 different random number generator 
seeds, and between 8 x 10 5 and 1.8 x 10 6 steps for acquisition were taken, after 10 4 steps for thermalization, on a 
system of size L = 64. Here a unit of time is considered to be one Monte Carlo step, that is L 2 single spin update 
trials. 

In Fig. U we show the results for T = 11.0,4.5,2.269,1.701,1.110, where temperatures are expressed in units of 
J/ks- We observe a two step decay also for very high temperatures. For all the temperatures except T = 11.0 we 
fit only the long time tail of the relaxation functions. For temperatures T > T p — 1.701 we could fit the long time 
regime a pure exponential, that is /3 = 1 within the errors in Eq. @. On the other hand for T = 1.110 < T p the 
long time behavior is not purely exponential, but can be fitted asymptotically with a stretched exponential. In Fig. 
H we show the values of /3(T) in function of T/T p . 



III. THE "g-BOND FRUSTRATED PERCOLATION" MODEL 

Using the Kasteleyn-Fortuin |2l[] and Coniglio-Klein |^2| cluster formalism for frustrated spin Hamiltonians |^3| , it 
is possible to show that the partition function of the model Hamiltonian (^]) is given by 

Z = Y^e^ n{c) q N{c \ (3) 
c 

where q = 2 is the multiplicity of the spins, (3 = l/fc^T, /i = kBT\n(e ql3J — 1), n(C) and N(C) are respectively 
the number of bonds and the number of clusters in the bond configuration C. The summation ^ c extends over 
all the bond configurations that do not contain a "frustrated loop" , that is a closed path of bonds which contains 
an odd number of antiferromagnetic interactions. Note that there is only one parameter in the model, namely the 
temperature T, ranging from to oo. The parameter fi, that can assume positive or negative values, plays the role of 
a chemical potential. 
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Varying q we obtain an entire class of models differing by the "multiplicity" of the spins, which we call the q- 
bond fully frustrated percolation model. More precisely, for a general value of q, the model can be obtained from a 
Hamiltonian |3l| 

U = -sjYy&iSiSi + 1)<W 3 - 2], (4) 

(y) 

in which every site carries two types of spin, namely an Ising spin and a Potts spin Ui = 1, . . . , s with s = q/2. For 
q = 1 the factor q N ( c ) disappears from Eq. (^), and we obtain a simpler model in which the bonds are randomly 
distributed under the conditions that the bond configurations do not contain a frustrated loop. For q — » we recover 
the tree percolation, in which all loops are forbidden, be they frustrated or not [3^ ]. 

When all the interactions are positive (i.e. ey = 1) the sum in (|J) contains all bond configurations without any 
restriction. In this case the partition function coincides with the partition function of the ferromagnetic g-state Potts 
model, which in the limit q = 1 gives the random bond percolation. 

From renormalization group and numerical results we expect that the model (|3|) exhibits two critical points p9| , p3| . 
The first at a temperature T p (q), corresponding to the percolation of the bonds on the lattice, in the same universality 
class of the ferromagnetic <7/2-state Potts model. The other at a lower temperature Tglq), in the same universality 
class of the fully frustrated Ising model. In the bidimensional case To(q) — (25-27 . 



IV. MONTE CARLO DYNAMICS 

A particular configuration of the model defined by (||) is determined by the state of each edge between two nearest- 
neighbor sites, that can be empty or occupied by a bond. The dynamics of the model is carried out in the following 
way: (i) choose at random a particular edge on the lattice; (ii) calculate the probability P of changing its state, that 
is of creating a bond if the edge is empty, and of destroying the bond if the edge is occupied; (iii) change the state of 
the edge with probability P. 

The point (ii) needs the knowledge of a non local property, namely if a bond placed on the chosen edge closes a 
loop or not, and if the loop is frustrated or not. This is accomplished in the following way: starting from the two 
sites at the extremes of the edge, visit the clusters of sites connected to them by a continuous path of bonds; if the 
clusters collide the bond closes a loop, otherwise it does not. Taking track of the number of antiferromagnetic bonds 
traversed visiting the cluster, one can determine also if the loop is frustrated or not. 

Note that at high temperatures clusters are small, and are visited in a few iterations, while at low temperatures 
density of bonds is high, and the clusters collide in a few iterations as well. On the other hand at the percolation 
transition clusters are very ramified, and one often must visit a large number of sites before the iteration is over. This 
makes the algorithm CPU consuming at the percolation transition, and prevents the simulation of very large systems. 

From Eq. (Q), the statistical weight of a bond configuration C is given by W(C) = e ^ n< - c \ N ^ if C does not 
contain a frustrated loop, and W(C) = if it does. Thus the transition probabilities for the principle of detailed 
balance must satisfy 

P{C -» C) = P{C -> C)e^ 5n q SN , (5) 

where Sn = n(C') — n(C) and SN = N{C) — N(C). Note that from Euler's equation SN = 5k — Sn, where 
Sn = ft(C') — k(C), and k(C) is the number of loops in configuration C, we can calculate SN knowing Sn and Sn. 
One can easily see that a possible choice for the transition probability P(C — > C) is given by 

„,„ _ / min(l,e /3 ^ 5 'y A ') if C is not frustrated, , , 

n ° ° > - { if C is frustrated. (bj 

The procedure described above in the points (i) — (iii) is called a "single update trial" . A Monte Carlo step consists 
in Q single update trials, where Q is the total number of edge on the lattice, namely, on the square bidimensional 



lattice, Q = 21? . In Sect. VI and VII, when we plot relaxation functions of the fully frustrated and locally frustrated 
bond percolation model, a unity of time is considered to be Q(p)~ 1 single update trials, or Monte Carlo steps, 

where (p) is the average density of bonds, ranging in the interval (0, 1). 
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V. STATIC PROPERTIES 



In this Section we analyse the percolation properties of the model defined by Eq. Q3J) , for q = 1 and q — 2, on a 
tridimensional square lattice, with fully frustrated interactions. We have used the histogram method for analyzing 
data H|. For each value of q, we have simulated the model for lattice sizes L = 32,48,64. For each size we have 
considered 10 temperatures around the percolation point, taking 10 4 steps for thermalization and between 3 x 10 5 
and 8 x 10 5 steps for acquisition of histograms. At every step we evaluate the following quantities: density of bonds 
p; existence of a spanning cluster P^; mean cluster size \. 

The quantity P^ assumes the value 1 if the bond configuration percolates, and if it does not. The mean cluster 
size is defined as |}5|] 

where n s is the number of clusters having size s on the lattice, and M — L 2 is the number of sites. The histogram 
method allows to evaluate the thermal averages of this quantities over an entire interval of temperature. The average 
of the quantity Poo is the probability of occurrence of a spanning cluster. 

Around the percolation temperature, the averaged quantities Poo(T) and x(P), for different values of the lattice 
size L, should obey the finite size scaling ]3(J 

P 00 (T)=F 00 [L 1 ^(T-T p )], (8a) 



X(T)=L^F X [LV"(T-T P )], (8b) 

where 7 and v are critical exponents of mean cluster size and connectivity length, P^ and F x are universal functions. 
Thus we can fit the values of v, 7 and T p so that plotting Pqo(T) and L-~t/ v x{T) in function of L l / V (T - T p ), the 
functions corresponding to different values of L collapse respectively on the universal master curves y = Foo[x\ and 
y = F x [x], Figures [| and || show the data collapse obtained for the q = 1 and q = 2 models, for lattice sizes 
L = 32, 48, 64. 

The values of the critical exponents extracted from the fit coincides, within the errors, with those of the ferromagnetic 
Potts model with spin multiplicity q/2 p^ |. Results are summarized in Tab. [|, while in Tab. [H| we report the critical 
exponents and transition temperature of the ferromagnetic Potts model. 



VI. THE RELAXATION FUNCTIONS OF THE "g-BOND FRUSTRATED PERCOLATION" MODEL 

In this Section we analyse the dynamical behavior of the model defined by (pi), simulated by the algorithm described 



in Sect. IV. For each temperature T and value of q, 32 different runs were made, varying the random number generator 
seed, on a system of size L = 64. We took between 10 3 and 10 4 steps for thermalization, and between 10 5 and 10 6 
steps for acquisition, calculating at each step the density of bonds p(t). The relaxation function of the density of 
bonds is defined as 

where 6p(t) = p(t) — (p). For each value of T and q, we averaged the 32 functions calculated, and evaluated the error 



as a standard deviation of the mean. As we mentioned in Sect. IV, we consider a unit of time to consist of Q{p) 
single update trials, where Q = 2L 2 is the number of edges on the lattice. 

In Fig. I the results for q = 1, T = 1.440, 1.067, 0.801, 0.625 are shown. Note that T p = 1.067 corresponds to the 
percolation transition of the model. For T > T p we fitted the calculated points with the function 

/(i)-/ exp-(t/r)' 3 . (10) 

The value of (3 extracted from the fit is equal to one within the error. Thus for T > T p the relaxation is purely 
exponential. For T < T p , we observe a two step decay, and only the long time regime of the relaxation functions could 
be fitted by Eq. ([h]). The value of (3 extracted is less than one, showing that stretched exponential relaxation has 
appeared for these temperatures. In Fig. [7| the values of 0(T) in function of the the ratio T/T p are shown, with least 
squares estimation errors. 
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In Fig. | the results for q = 2, T = 2.269, 1.701, 1.440, 1.110 are shown. Temperature T p = 1.701 corresponds to the 
percolation transition. The fits were made in the same way described for q = l, and the values of (3(T) extracted are 
shown in Fig. ^. Also in this case = 1 within the error for T >T p , and (3 < 1 for T < T p . Note that the q = 2 fully 
frustrated bond percolation model can be mapped exactly onto the fully frustrated Ising spin model, as we showed 



in Sect. III. So it is interesting to compare the relaxation functions of this model to those of the corresponding fully 



frustrated Ising model, simulated by standard spin flip techniques. 



VII. THE LOCALLY FRUSTRATED BOND PERCOLATION 



In the fully frustrated g-bond percolation model, the configurations of bonds which contain at least one frustrated 
loop have zero weight. The size of frustrated loops has no upper limit. To study systematically the effect of frustration, 
we consider now a modified version of the model, in which only loops up to some specified length are considered 
frustrated, while longer ones are permitted. The partition function of this model is given by 

c 

Here the parameters j3 and fi have the same meaning as in Eq. (^), n(C) is the number of bonds in the configuration 
C, and the sum 53c ^ s extended over the bond configurations that do not contain frustrated loops of length 
A < Ao • Thus the model switches continuously from the random bond percolation (Ao = 0) , and the g-bond frustrated 
percolation with q — 1 (Ao = oo). 

We have studied the critical properties of the model, for Ao = 4, and its dynamical behavior above and below the 
percolation transition. The critical exponents extracted from the finite size scaling data collapse (see Fig. |Io| ) coincides 
within the errors with those of the random bond percolation, as shown in Tab. [j]. The percolation temperature is 
T p = 1.277, intermediate between that of the random bond percolation, T p = 1.443, and that of the q = 1 bond 
frustrated percolation, T p — 1.067. 

The dynamics, on the other hand, can feel the local constraint constituted by the frustration, as the autocorrelation 
functions do not decay as a single exponential when the temperature is lowered below the percolation threshold. 
However, the long time regime of the relaxation functions could be fitted with an exponential for all the temperatures 
considered. 



In Fig. 11 we show the bond density autocorrelation functions, evaluated on a lattice L = 64, for temperatures 
T = 1.277,0.911,0.625. The function calculated for T — 1.701 is not shown because it overlaps with the function 
calculated for T = 1.277. Solid curves show fits made by a pure exponential. Averages were taken over 32 different 
runs, each one taking 10 4 steps for thermalization, and between 2 x 10 5 and 3 x 10 5 steps for acquisition. 



VIII. CONCLUSIONS 



We have studied the fully frustrated Ising model and the fully frustrated percolation model. The dynamics of the 
models was analysed in detail, in connection with the problem of the onset of stretched exponentials in frustrated 
systems, like glasses and spin glasses. 

Due to absence of disorder in these models, the arguments suggested by Randeria et al, which predict a non- 
exponential relaxation below the critical temperature T c of the corresponding ferromagnetic model, does not apply. 
In fact our results show no sign of complex dynamical behavior at T c . We find instead an exponential relaxation 
above the percolation temperature T > T p , while for T < T p the long time tail of the relaxation functions can be 
fitted with a stretched exponential. So we conclude that at least in the models without disorder the appearance of 
complex dynamics is related to a percolation transition. In systems like spin glasses T p and T c are very close, and it 
is difficult to distinguish numerically where the onset of non exponential relaxation occurs. 

We also find that frustration plays an important role in the critical properties at the percolation threshold T p . For 
example, for q — 1 the critical behavior is in the same universality of the ferromagnetic q = 1/2-state Potts model, 
contrary the unfrustrated case, which corresponds to random bond percolation and is in the same universality class 
of the q = 1 ferromagnetic Potts model. 

We have also considered a model, the locally frustrated bond percolation, in which only loops up to length four 
are considered frustrated. The model shows the same critical properties as the random bond percolation, showing 
that the frustration is "too local" to change the universality class. Similarly, the relaxation functions in the long time 
regime can always be fitted with an exponential, showing that the frustration constraint is not enough to give rise to 
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stretched exponential relaxation. More careful study of this model, possibly varying the "range" of the frustration 
between the size of the single plaquette, and that of the whole system, may shed more light on the role played by the 
frustration in the dynamics of complex systems. 
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FIG. 1. Distribution of interactions for the fully frustrated model. Straight lines and wavy lines correspond, respectively, to 
eij = 1 and eij = — 1. 



FIG. 2. Relaxation functions f(t) of energy in function of time t for the fully frustrated Ising model, with spin flip dynamics, 
lattice size L = 64, for temperatures (from left to right) T=ll. 0,4.5, 2. 269, 1.701, 1.110. 
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FIG. 3. Stretching exponents j3(T) in function of T/T p , the ratio of temperature over percolation temperature, for the fully 
frustrated Ising model, with spin flip dynamics, lattice size L = 64. 



FIG. 4. Finite size scaling of (a) P tx ,(T) and (b) x(T), for the q = 1 model, and for lattice sizes L — 32, 48, 64. 



FIG. 5. Finite size scaling of (a) Poo(T) and (b) x(T), for the q = 2 model, and for lattice sizes L — 32, 48, 64. 

FIG. 6. Relaxation functions f(t) of bond density in function of time t for q = 1, lattice size L — 64, for temperatures (from 
left to right) T=1.440,1.067,0.801,0.625. 

FIG. 7. Stretching exponents f3(T) in function of T/Tp, the ratio of temperature over percolation temperature, for the q = 1 
fully frustrated bond percolation model, lattice size L = 64 

FIG. 8. Relaxation functions f(t) of bond density in function of time t for q — 2, lattice size L = 64, for temperatures (from 
left to right) T=2.269,1.701, 1.440, 1.110. 

FIG. 9. Stretching exponents (3(T) in function of T/T p , the ratio of temperature over percolation temperature, for the q = 2 
fully frustrated bond percolation model, lattice size L = 64 



FIG. 10. Finite size scaling of (a) Poo(T) and (b) x(T), for the local frustrated model, and for lattice sizes L — 32, 48, 64. 

FIG. 11. Relaxation functions fit) of bond density in function of time t for the local frustrated model, lattice size L = 64, 
for temperatures (from left to right) T=l. 277,0. 911,0. 625. 



model 


l/i/ 


7 


T P 


q = l 


0.56 ±0.02 


3.22 ±0.07 


1.067 ±0.001 


q = 2 


0.75 ±0.03 


2.34 ±0.06 


1.701 ±0.001 


local (A =4) 


0.75 ±0.03 


2.33 ±0.04 


1.277 ±0.001 



TABLE I. Critical exponents 1/v and 7, and percolation transition temperature T p of the fully frustrated q-bond percolation 
model with q = 1,2, and of the locally frustrated bond percolation model with A = 4. 



model 


1/u 


7 


T c 


9=1/2 


0.56 


3.27 


1.233 


q = l 


0.75 


2.39 


1.443 


q = 2 


1 


1.75 


2.269 



TABLE II. Critical exponents 1/u and 7, and critical temperature T c of the ferromagnetic Potts model, with multiplicity of 
spins q = 1/2, 1, 2. 
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